arXiv:1507.03871v2 [astro-ph.CO] 30 Oct 2015 


APCTP-Pre2015-019, RUNHETC-2015-06 


WIMP isocurvature perturbation and small scale structure 

Ki-Young Choi,i’Q Jinn-Ouk Gong,^- 3-Q and Chang Sub Shinb[^ 

^Korea Astronomy and Space Science Institute, Daejeon 305-348, Korea 
^Asia Pacific Center for Theoretical Physics, Pohang, 790-784, Korea 
^Department of Physics, Postech, Pohang 790-784, Korea 
^Department of Physics and Astronomy, Rutgers University, Piscataway NJ 08854, USA 

The adiabatic perturbation of dark matter is damped during the kinetic decoupling due to the 
collision with relativistic component on sub-horizon scales. However the isocurvature part is free 
from damping and could be large enough to make a substantial contribution to the formation of 
small scale structure. We explicitly study the weakly interacting massive particles as dark matter 
with an early mater dominated period before radiation domination and show that the isocurvature 
perturbation is generated during the phase transition and leaves imprint in the observable signatures 
for small scale structure. 

PACS numbers: 95.35.+d, 14.80.Ly, 98.80.Cq 


Introduction. The formation of large scale structure 
is consistent with non-relativistic dark matter (DM) in¬ 
dependent of its nature. Small scale structure, however, 
depends on the microphysics of DM and the correspond¬ 
ing evolution in the early universe W- For weakly in¬ 
teracting massive particles (WIMPs), the kinetic decou¬ 
pling is a crucial stage to determine the size of smallest 
object during the process of kinetic decoupling col- 

lisional damping smears out the inhomogeneities below 
the corresponding damping scale. After kinetic decou¬ 
pling WIMPs can move freely and this leads to additional 
damping below the free streaming scale. For neutralino 
DM, the kinetic decoupling scale is set when the temper¬ 
ature is 10 MeV - 1 GeV for the mass between 100 GeV 
and TeV [7]. 

In radiation dominated era (RD), while the “adiabatic” 
component of DM perturbation on sub-horizon scales ex¬ 
periences oscillations followed by collisional damping [5] , 
the isocurvature perturbation between DM and radia¬ 
tion, 


=6^-hr, (1) 

\ Pm Pr J 4 

remains constant without damped oscillations 0 HU]. 
This property was used to explain large scale structure 
with baryon isocurvature perturbation 0 , which is ruled 
out now by the adiabatic constraint from the cosmic mi¬ 
crowave background (GMB) [TT]. However, large isocur¬ 
vature perturbation on small scales is not constrained by 
the GMB observations and can give observable signatures 
in small scale structure. 

In this article, we show how large isocurvature pertur¬ 
bation of WIMPs can be generated for scales that enter 
the horizon before the kinetic decoupling. If S' = 0 at the 
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onset of RD, it remains so during kinetic equilibrium. In¬ 
stead, if an early matter dominated era precedes RD, a 
sizable amount of S can be generated. We note that this 
isocurvature perturbation will not be damped even if the 
kinetic decoupling happens after the transition to RD. 

Dark matter in non-thermal background. In the early 
universe, it happens often that the energy density of the 
universe is dominated by a non-relativistic matter which 
subsequently decays into relativistic particles. This non- 
relativistic matter includes a coherently oscillating scalar 
field like an inflaton, or massive fields which decay very 
late, such as curvaton, moduli and so on. As an illustra¬ 
tion, we consider this dominating non-relativistic mat¬ 
ter as a scalar (p with a decay rate F,^. Accordingly, 
we call the epoch during which cp dominates the energy 
density as the scalar dominated era (SD). In the back¬ 
ground, then there are three species of fluid: </>, radiation 
and DM. Their evolutions are governed by the continuity 
equations. 
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where M is the mass of the DM particle, fm is the frac¬ 
tion of the decay of p into DM, (aav) is the thermal 
averaged annihilation cross section of DM and pff} « 
M'^(27rM/r)“3/2 exp(—M/T) is the energy density of 
DM in thermal equilibrium. Here radiation is the rela¬ 
tivistic particles thermalized quickly when produced from 
the decay of p, and thus the temperature T is prop¬ 
erly defined by its energy density pr = with 

< 7 * being the effective degrees of freedom of the rela¬ 
tivistic particles in thermal equilibrium. The reheat¬ 
ing temperature is then approximately given by T^eh ~ 
(7r^5*/9O)“3/4\/^pjF0. For successful big bang nucle¬ 
osynthesis, we require that T^eh must be larger than 
C>(MeV) [II]. 
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FIG. 1: The evolution of the energy densities of the scalar 
(black), radiation (red) and DM (blue) respectively with re¬ 
spect to the initial total energy density. Blue dashed line is the 
equilibrium energy density of WIMP, and green dashed lines 
denote the asymptotic behavior of radiation energy density. 
DM freezes out at ajai ~ 20 and RD starts from ajat ~ 300. 

While radiation is produced directly from the decay of 
(p, DM can be produced in several different ways [l3] . For 
simplicity, we assume that DM is produced only from ra¬ 
diation by scatterings and set fm = 0. Even in this case, 
a sizable amount of DM can be produced from thermal 


plasma. If the interaction of DM with plasma is large 
enough, they could be in thermal equilibrium. WIMP 
is one such example, which is intimately coupled to the 
relativistic plasma and decoupled when T/M ^ 1/20, de¬ 
pending on the annihilation cross section {aav) |14j . The 
freeze-out may happen during SD or RD after the scalar 
decay. For the latter case, there will be no difference from 
the thermal WIMP in the standard scenario. Therefore, 
in our study, we will focus on the case that WIMPs are 
decoupled during SD. 

In Figure we show the evolution of the background 
energy densities of (p, radiation and DM by solving (H- 
During SD, pr scales as pr oc a due to the con¬ 
tinuous production from the scalar decay and thus the ef¬ 
fective equation of state during SD is —1/2. DM is frozen 
during SD, and its energy density decreases simply pro¬ 
portional to a~^ after then. However the interactions by 
collisions continue until RD. 

Evolution of perturbations. Now we consider the evo¬ 
lution of perturbations. For this, we use the Newtonian 
gauge with the metric 

ds^ = —(1 -I- 2$)dt^ -I- a^(l — 2^)5ijdx'‘dx^ . (5) 

The perturbation equations can be derived from the 
Boltzmann equation for each component {a = (p, r and 
m) and they are given by [ails] 
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where 6a = V ■ Va = divl, is the velocity divergence field, 
= Wm = 0 and Wr = 1/3. At leading order of T/M, 
the energy-momentum transfer functions Qa and SQa 
can be calculated from the Boltzmann equation as 
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where we have put fm = 0. In the above equations, 
we have included the elastic scattering cross section be¬ 
tween radiation and DM ae which keeps DM and radi¬ 
ation in kinetic equilibrium until they decouple at 
set by Ce{crev)pr/M\T=Tt,i = 17(Fkd), with Ce = 0(1) be- 
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FIG. 2: The evolution of the density contrast of the radiation (red), DM (blue) and the isocurvature perturbation (brown) with 
respect to the initial gravitational potential for < k~^ {k = 0.1 fckd, left), < k~^ < k^J^ {k = 5fckd = O.Sfcreh, middle) 
and < k~^ < k~l (k — SOfckd = 5 fereh = O.Sfcu, right). We have set M = 5TeV, T-eh = 0.1 GeV and Tkd = 0.01 GeV. 


ing dependent on model due to the different momentum 
dependence of CTg. 

The 00 component of the perturbed Einstein equation 
governs the evolution of the metric perturbations, 


— lb — 3H flit + 77$^ — ^ + PrSr + PmSm) ■ 

a ^ / .^777/pi 

(17) 

In the absence of the anisotropic tensor, we can set $ = d' 
which then closes the above set of equations. This is 
possible since (j) and radiation which dominate the en¬ 
ergy density are isotropic in our setup. Note that the 
effects of the anisotropic shear and non-vanishing sound 
speed of DM, Cg ^ y/T/M, can be important after ki¬ 
netic decoupling for scales smaller than the free stream¬ 
ing length In m, it is shown that when the free 

streaming length is much shorter than the scale that 
enter the horizon at the moment of kinetic decoupling, we 
can take an approximation that solving the Boltzmann 
equations first in perfect fluid limit while maintaining 
the elastic scattering, and then multiplying the solution 
by the Gaussian suppression term. Actually this limit is 
also physically interesting, because two different damping 
scales can be more clearly distinguished. 

In this article, we consider the hierarchies among scales 
as k^^ < , where is the scale that en¬ 

ters the horizon at T = Tj-eh- This means that the free 
streaming scale enters the horizon during SD and that 
kinetic decoupling occurs during RD. The large hierar¬ 
chy between and k^^ can be obtained when M is big 
enough while the elastic scattering is mediated by a field 
much lighter than DM. In this case, the freeze-out abun¬ 
dance also could be large, but the subsequent dilution by 
entropy injection from the scalar decay can provide the 
correct amount of the present DM density [T71IT5] . For 


WIMP, we find [T^ 



where g^s is the effective number of light species for en¬ 
tropy and Teq = ©(eV) is the temperature at matter- 
radiation equality. 

In Figure we show the evolution of perturbations 
on three different scales. During SD, the perturbations 
are adiabatic on super-horizon scales since both radiation 
and DM are produced from a single source (j), which set 
the initial values of perturbations as 5^{ai) = 25r{ai) = 
—and 5m{a.i) ~ M5r{ai)/{ATi), with Ti being deter¬ 
mined from pr{ai). During the transition from SD to RD, 
$ rescales from to 10<I>i/9 on super-horizon scales and 
accordingly Sr changes from —to —2(10/9)$^. Mean¬ 
while, at early times when DM is in thermal (chemical) 
equilibrium, Sm oc and is reduced to —5<I>i/3 during 
RD which follows the adiabatic condition 6m = 3(5r/4. 

While for modes which enter the horizon after ki¬ 
netic decoupling (k^^ < k~^), Sr oscillates and Sm 
grows logarithmically as shown in the left panel of Fig¬ 
ure!^ for the modes which enter before kinetic decoupling 
(^reh < < k^^) Sm oscillates together with Sr and is 

damped, which is known as collisional damping. The 
non-vanishing sub-horizon entropy perturbation appears 
due to the damping of Sm as shown in the middle panel 
of Figure 

An interesting feature happens for the modes which 
enter the horizon during SD but after the free stream¬ 
ing scale enters (k^^ < k~^ < k~f^) as in the right 
panel of Figure During the transition from SD to 
RD, Sm does not follow Sr, and the isocurvature per¬ 
turbation is generated. In this period, DM is no longer 
produced after chemical freeze-out and the number den¬ 
sity is frozen while radiation is still being produced from 
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(j). The continuous entropy injection becomes the source 
of the isocurvature perturbation between DM and radi¬ 
ation. This perturbation still persists even after kinetic 


decoupling. Before calculating its analytic expression, we 
explicitly show why it is not damped, from the solution 
for Sm during RD [IB] . 
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where and 4>k are fc-dependent constants while Ak{t) 
and Bk(t) vary in time. Their time dependence is deter¬ 
mined by the elastic scattering term as 
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The values of A^, (pk, ^fc(ireh) and i?fc(treh) are given at 
the onset of RD, and for adiabatic modes they are 


Afc = —10<i)j, (j)i.=0, Afe(treh) = — 104)i, 

Rf (treh) = -10$. [lE - I 


(23) 


where je ~ 0.577 is the Euler-Mascheroni constant. 
Then on super-horizon scales k aH we can recover 
—5$i/3 during RD. For the modes which enters during 
RD (ic“(j < fc“^), the solution is [16] 
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for (aev) oc T^+", which clearly shows the damping for 
k~^ ^ k^^ due to the collision with radiation. 


Here it is important to note that in (22) only Bk ap¬ 
pears. The additional constant term to the adiabatic 
one is not damped away even in the kinetic equilibrium 
/ decoupling periods. As a result, for k~^ k^^, Sm 

is dominated by the isocurvature perturbation: Bk = 


Br 


^ad 


Generation of isocurvature perturbation. For the 
modes that enter the horizon during SD after chemical 
decoupling of DM, S^ grows linearly, 


S^a) = -2$. - 
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and then logarithmically during RD. Meanwhile, Sr grows 
during SD, since radiation is continuously produced from 
the decay of 4>. However, after the transition from SD 
to RD, this enhancement is lost and Sr oscillates with 
heavily suppressed amplitude [T|. 

During kinetic equilibrium, DM is tightly coupled to 
radiation, so that 9m ~ Ignoring the effect of DM 
annihilation the relevant equations for Sm and Sr are, 
from 


a 

Sr^-^- + ^{S^-Sr) , 

6 a pr 
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where we have neglected 0(1) contribution. From SD to 
the transition period, both Sr and $ are sub-dominant 
compared to S^, and pr ~ 2Trj,prj,/5H. Then the isocur¬ 
vature perturbation is 


*S'(^reh) ~ ^ 


dt 


^,pp<pS^ 




k 
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As can be read from (27), unlike Sm, Sr is sourced by both 


Or and S^ because there is steady production of radiation 
from (p. The corresponding isocurvature part becomes 

. . 

While the isocurvature perturbation can avoid the 
damping due to the collision, the diffusion by the free 
streaming still exist. Considering the damping effect due 
to free streaming, as discussed before we may add a Gaus¬ 
sian suppression factor to Sm as 


exp - 


2kl 




—y 

^reh / 


(29) 


where the free streaming scale is estimated as (201. 


Based on these results, it is straightforward to calculate 
the evolution of perturbation during the subsequent mat¬ 
ter dominated era, the transfer function and the mass 
function. In Figure we show Sm at later stages at 
a = 100Okd and lOoeq. 

Implications. At the heart of our finding is that de¬ 
spite of the damping of the conventional adiabatic per¬ 
turbation of DM due to a large elastic scattering rate 
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FIG. 3: Density contrast of DM with M = 5TeV, Treh = 
0.1 GeV and Tkd = 0.01 GeV. 

between DM and the standard model particles, the DM 
isocurvature perturbation survives the collisional damp¬ 
ing until kinetic decoupling. This unsuppressed pertur¬ 
bation on small scale can give rise to a large number of 
DM clumps, such as compact mini haloes [201 ED- 
Since DM can annihilate efficiently in the clumps, 
these haloes can serve as the sources of highly luminous 
gamma rays which can be well observed with the ongo¬ 
ing or future gamma-ray telescope like Fermi-LAT [22] 
or Cerenkov Telescope Array [23]. They can be also the 


sources of neutrinos [23] , detectable by IceCube [23] . Fur¬ 
thermore they can leave an imprint in the CMB by chang¬ 
ing the reionization history of the Universe |26j . produce 
a microlensing light curve [27] [2H], or change the direct 
detection rate [23]. 

The requisite for a large enough DM isocurvature per¬ 
turbation is a sufficient hierarchy between DM freeze-out 
and reheating to have a long enough early MD. This is 
easily realized with a low reheating temperature, which 
happens ubiquitously in many theoretical models when 
the heavy particle dominates and decays in the early Uni¬ 
verse. Those models include the neutralino DM in the 
low reheating temperature [SO] El] and the scenario of 
decaying heavy particle such as moduli, gravitino [32] j or 
axino [18] . Considering both astrophysical and cosmo¬ 
logical observations and DM theories should give more 
information about the early history of the universe be¬ 
fore BBN and the properties of DM. 
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